Quantum homodyne tomography of a two-photon Fock state 
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We present a continuous-variable experimental analysis of a two-photon Fock state of free-propagating 
light. This state is obtained from a pulsed non-degenerate parametric amplifier, which produces two intensity- 
correlated twin beams. Counting two photons in one beam projects the other beam in the desired two-photon 
Fock state, which is analyzed by using a pulsed homodyne detection. The Wigner function of the measured 
state is clearly negative. We developed a detailed analytic model which allows a fast and efficient analysis of 
the experimental results. 
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Quantum properties of light beams can be described in 
terms of amplitude and phase or, in Cartesian coordinates, in 
terms of the "quadrature components" of the quantized elec- 
tric field, associated with non-commuting operators x and p. 
The corresponding observables, often called "quantum con- 
tinuous variables", are analogous to the position and the mo- 
mentum of a particle, and from Heisenberg's inequalities they 
cannot be determined simultaneously with an infinite preci- 
sion. As a consequence, one cannot define a proper phase- 
space distribution I{{x,p) for the electric field, but rather a 
quasidistribution W{x^p) called the Wigner function. This 
function can be reconstructed by quantum homodyne tomog- 
raphy (U, which consists in measuring several quadratures 
xg = X COS 6 + p sin 6 with a homodyne detection, and 
applying an inverse Radon transform. 

The most conspicuous property of the Wigner function is 
that it may take negative values for specific quantum states, as 
a signature of their non-classical nature. This is the case for 
Fock states, which contain a well-defined number of photons. 
Such states can be generated by using "twin" beams, which 
are produced by optical parametric amplification, and which 
contain perfectly correlated numbers of photons. Counting n 
photons in one mode projects the other mode in a n-photon 
Fock state, which can then be analyzed using homodyne to- 
mography. This was recently demonstrated for n = 1 J^IH. 
However, up to now this method could not be applied for 
higher photon numbers, since the probability to generate si- 
multaneously more than one photon pair was extremely low. 

In this Letter we present a detailed analysis of a free- 
propagating light pulse prepared in a two-photon Fock state 
{n = 2). The measured Wigner function presents a complex 
structure and takes negative values. In addition to standard 
methods, we will also present a novel analytic model of the 
experiment, allowing an in-depth physical interpretation of the 
experimental results. 

Our experimental setup is presented on Fig. ^ A pulsed 
Ti-Sapphire laser produces 180-femtosecond nearly Fourier- 
limited pulses with an energy of 40 nJ and a 800 kHz repe- 
tition rate |4]. The high pulse peak power allows us to in- 
crease thepair production rate beyond what was available pre- 
viously The 850 nm pulses are frequency-doubled 
[second harmonic generation (SHG)] by a single pass in a 



100 /im thick non-critic ally phase-matched potassium nio- 
bate (KNbOa) crystal. The frequency-doubled beam pumps 
an identical crystal used as an optical parametric amplifier 
(OPA), generating a two-mode squeezed state 15]. To align 
the setup, a probe beam is injected in the OPA with an angle 
of 5° to the pump direction. It allows to measure a classical 
phase-independent gain g = 1.07. The homodyne detection is 
aligned on the idler beam, whereas the signal beam, after spa- 
tial and spectral filtering, is split between two avalanche pho- 
todiodes (APD) operating in a photon-counting regime. The 
detection of a coincidence by the APDs means that at least 
two photon pairs were created in the OPA by the same pulse. 
Since the gain g is still relatively low the probability to create 
more than two pairs is small in this case. Therefore, a coin- 
cidence detected by the APDs conditionally prepares a two- 
photon state in the idler beam. Single-photon states are condi- 
tioned by single APD events. The prepared states are analyzed 
by a homodyne detection operating in a time-resolved regime. 
It samoles each individual Dulse. measurine one Quadrature 




FIG. 1 : Experimental setup, and Wigner function of the two-photon 
state propagating in the experiment (corrected for homodyne detec- 
tion losses, see text). 
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TABLE I: Critical values of the Wigner functions corresponding to 
the measured uncorrected data (Raw, obtained from the Radon trans- 
form), to the state corrected for homodyne detection losses (Cor- 
rected, obtained from the MaxLik method) and to the ideal state 
(Ideal). 





2 photons 


1 photon 




min{W2) W2{Q) 


min{Wi) = Wi{0) 


Raw 


-0.009 ±0.003 0.012 ±0.003 


-0.052 ± 0.003 


Corrected 


-0.034 ± 0.003 0.062 ± 0.003 


-0.123 ±0.003 


Ideal 


-0.13 0.32 


-0.32 



FIG. 2: Experimental quadrature measurements, and quadratures re- 
constructed using our model (see text) 
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FIG. 3: Wigner functions of the measured uncorrected states, recon- 
structed by a standard Radon transform, compared to those obtained 
using the model described in the text 



Xq in phase with the local oscillator. 

In previous n — 1 state reconstruction experiments 
it was generally admitted that the generated states are phase- 
independent. In our case, the production rate of single photons 
is very high, and we can record the full n — 1 quadrature 
distribution in less than a second, during which phase drifts 
are negligible. Therefore we did check experimentally that 
both the unconditional (thermal) and singly-conditional {n = 
1) probability distributions do not depend on 9. Then it is 
quite reasonable to assume that this is also the case for the n = 
2 state, as it was done for the n = 1 state in older experiments. 

In a 2-hour experimental run we acquired 105.000 ho- 
modyne data points conditioned on two-photon coincidences 
(40 seconds were enough to acquire 180.000 single-photon 
events). Dividing the data into 64-bin histograms, we obtained 
the quadrature distributions presented on Fig. |2] With a nu- 
merical Radon transform, we reconstructed the Wigner func- 
tions associated with the measured states (see Fig. 13, both 
clearly negative. Their minima and their values at the origin 
are presented in Table U To determine the Wigner functions 
of the generated states, presented on Fig. |3 we correct for 
the homodyne detection losses by using a standard maximal- 
likelihood (MaxLik) algorithm IiIIsIIj taking into account an 
independently measured homodyne efficiency = 80%. 

The negativity of the Wigner function can be rapidly lost 
with experimental imperfections. Above all, we must ensure 
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FIG. 4: Experimental Wigner functions corrected for losses in the ho- 
modyne detection, reconstructed by a standard maximal-likelihood 
(MaxLike) method, compared to those obtained using the model de- 
scribed in the text 



that the prepared state belongs to the mode analyzed by the 
homodyne detection. This modal overlap ^ is decreased by 
the imperfections of the filtering system, by the APD dark 
counts, and by the limited spectral and spatial qualities of the 
optical beams. As a result, we may consider that the state is 
prepared in the right mode with a probability ^, and in an or- 
thogonal mode with a probability 1 — ^. A second source of 
decoherence is excess noise in the OPA, producing uncorre- 
cted photons. The actual OPA can be represented by an ideal 
non-degenerate amplifier with a gain g = cosh^(r), produc- 
ing a pure two-mode squeezed state, followed by two phase- 
independent amplifiers on signal and idler beams, each one 
with a gain h = cosh^(7r), where 7 is the ratio between the 
undesired and the desired amplification efficiencies (ideally 
7 = 0). Finally, the homodyne detection presents a finite 
efficiency 77 and an excess noise e. From the measured op- 
tical transmission rjt — 97%, quantum detection efficiency 
r]q = 97.5% and mode-matching efficiency r],n = 92%, we 
estimate 77 = rjtrjqT]^ = 80%. Since 7/ and e are not in- 
volved in the preparation but only in the analysis of the state, 
we can correct for their effects in order to determine the actual 
Wigner function of the generated state. The overall efficiency 
/i of the APD detection channel, although rather low (6%), is 
not a limitation in this experiment (see Appendix). 

In order to obtain a more physical analysis of our data, we 
have constructed a complete - but nevertheless simple - an- 
alytic model of the experiment (see Appendix). Apart from 
predicting the performance of the setup, it allows to extract 
much more information from the experimental data than the 
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numerical methods presented above, although it is, of course, 
less general. It uses a generic parameterized expression of the 
Wigner function, derived in the Appendix, which accounts for 
all the experimental defects : 



W2{x,p) 

where 



(1- 5^ + 2(1 -5)'-^ 



cr^ = 2Tj{hg - 1) + 1 + e 

S = 2^rjh^g{g-l)/[a\hg^l)] 



The associated quadrature distribution is described by 
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For the one-photon case, the same method leads to 
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The density matrices of these states are diagonal in the Fock 
basis, the non-zero coefficients given by : 

'^IpM = , .i^, [si - 2n{n + l)S'a'] (7) 



(^2 ^ i)n+3 

(n|pi|n> = 25„(r72 - 1)"- 7(^2 + 



(n\po\n) 



2(a' 



1)7(^2 



1 



(8) 
(9) 



where 5„ — (T'''(l—(5)+(T^(5(l+2ri) — l, and po corresponds 
to the thermal unconditioned state (obtained by taking 5 — 
in any of the above equations). 

These states are completely described by the two same pa- 
rameters (T^ and 5. Here is simply the variance of the non- 
conditioned gaussian thermal state. The non-classicality of 
the conditioned states is determined by 6, which varies be- 
tween for a non-conditioned state and 2 for the ideal case. 
When (5 > 1, both Wi and W2 become negative, and a cen- 
tral peak appears on W2. These parameters, very useful to 
optimize the experiment, can be directly extracted from the 
second and fourth moments of the measured distributions : 



1 photon 



2 photons 



{x^)i = a^l + 5)/2 {x^)2 - a^l + 2S)/2 

= 3a4(l + 2,5)/4 {x'^)2 = + 46 + S^)/A 

We used one-photon conditioning during the optimization, 
so that (7^ and 6 could be determined in a few seconds, 300 
times faster than in the two-photon case. The two-photon 
state, described in principle by the same parameters, was "au- 
tomatically" optimized in this process. We found that the val- 
ues deduced from single and two-photon state tomographies 
are exactly the same for a'^, and differ by less than two per- 
cent for S. 
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FIG. 5: Main density matrix coefficients of the states conditioned on 
0, 1 and 2 photons (corrected for detection losses). 



In addition, the quadratures reconstructed using the pa- 
rameters (T^ and S extracted from raw data are in excellent 
agreement with the measurements (see Fig. |2}, and the re- 
constructed Wigner functions of the measured states are very 
close to those obtained by the Radon transform (Fig.|3}- Equa- 
tions|2]and|3]also allow to determine the modal overlap ^ and 
the excess gain parameter 7. The obtained values (^ — 0.9 
and 7 = 0.4) are fully compatible with experimental evalua- 
tions, which are difficult to do but were carried out by using 
independent classical amplification and photon counting tech- 
niques. 

Since the results obtained with this method appear to be 
completely consistent, both within themselves and with inde- 
pendant measurements, we can assume that the Wigner func- 
tion of the generated state, which we would measure with an 
ideal homodyne detection, can be simply calculated by taking 
r] = \ and e = in our expressions, keeping all other parame- 
ters unchanged. The obtained results are again in good agree- 
ment with those provided by the maximal-likelihood method, 
as shown on Fig. 0] The main density matrix coefficients 
of the generated states are represented on Fig. |5] This gives 
confidence that our method provides a very fast and reliable 
way to interpret the experimental data, which is more "con- 
strained" than the Radon transform, but also much closer to 
the physics of the experiment. 

The present experimental and theoretical results demon- 
strate simple techniques to generate and analyze sophisti- 
cated non-classical states of propagating light fields, which 
have been considered almost out of experimental reach during 
many years. Similar methods can be used to create photon- 
subtracted entangled states with two-mode negative Wigner 
functions, which should improve the fidelity in teleportation 
experiments I9l \\(X |l III, and allow to implement loophole- 
free Bell tests pizUlSl]. The avenue of manipulating negative 
Wigner functions now seems clearly open for quantum com- 
munications. 



APPENDIX 

The model for the experiment is represented on Fig. |6l The 
OPA produces a two-mode noisy squeezed state with a density 
matrix psqz associated with a Wigner function 



Wsqz{xi,Pl,X2,P2) 



(10) 
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FIG. 6: Modelling of the experiment. 
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where s 



is the two-mode variance squeezing factor as- 



sociated with a gain g — cosh^(r), and h — cosh^(7r) is the 
excess gain. The mode 1 is directed towards the homodyne 
detection, whereas the mode 2 is sent into the conditioning 
channel. The homodyne losses can be represented by mix- 
ing the mode 1 with vacuum on a beam splitter (BS) with a 
transmission T = rj. Since we are only interested in the trans- 
mitted mode H, we trace over the reflected mode to obtain the 
resulting density matrix. The same holds for the APD losses, 
with a transmission T = ji. 

The resulting Wigner function Wioss is calculated by con- 
volution of Wsqz with the Wigner functions W^ac of two vac- 
uum modes, using Wyac{x,p) — exp(— — p^)/7r. Then, 
the mode transmitted through the APD channel is mixed with 
another vacuum mode on a 50/50 beamsplitter, producing a 
density matrix pmix involving three modes H, A and B, and 
associated with the Wigner function : 



W„ 



= Wlo.s ( 



xh,Ph, 



xa+xb PA+pe 



72 



72 



PA~PE 



The modes A and B are detected by the APDs A and B, which 
realize respectively the projective measurements Ha.b = 
Id — |0a.s)(0a,_b| with a probability ^ ("matched clicks"), 
and Ho = Id with a probability 1 — ^ ("unmatched clicks"). 
The density matrix becomes 

Pcond = N2i^'ilA^BPmix^A^B + (1 " O'^ Prmx 

+Ni£,{l - £,){UAPmix^A + ^BPnUx^s) 

where Ni = l/Tr{UAPrmx) = l/Tr(nBPrmx) and N2 - 
l/Tr{IlA^BPmix)- Finally, the density matrix of the mea- 
sured two-photon state is obtained by tracing out the two APD 
modes A and B : 

P2 = Tta.B Pcond 

= [N2e + 2iViC(l - C) + (1 - 0'] Tta^B Pra^x 



- {f^2e + ^1^(1 - 0] TrB{^A\p^vx\^A) 

- + ^1^(1 " 0] rr^(OB|p™„|OB) 

^N2i^l^^A^B\p^^x\^A^B) (H) 

The associated Wigner function can be calculated using 

TrxW^mx = J WrmxdxxdpK 

{OK\W,nix\OK) = StT J W^i^WyacdXKdpK 

where K — A,B. As expected, it has no definite phase and 
depends only on = + p\j. It has the form 



W. = 



^ /3e ^ (l-a + /3)e" 



(12) 



where a, (3 and (Xj are functions of the parameters above. This 
linear combination of gaussian functions looks quite simple, 
but a and (3 diverge when the OPA gain or the APD efficiency 
are small, which is our case. This leads to numerical instabil- 
ities when this expression is used for data analysis. To avoid 
this problem one can simply take the Umit p, ^ in eq. il2\ . 
obtaining eq. Q quoted in the main text above. In our range 
of parameters, these two equations are numerically indistin- 
guishable. 
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